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ABSTRACT 

We report evidence for an anti-correlation between spin temperature T s and metallicity [Z/H], detected at 
3.6cr significance in a sample of 26 damped Lyman-a absorbers (DLAs) at redshifts 0.09 < z < 3.45. The 
anti-correlation is detected at 3<7 significance in a sub-sample of 20 DLAs with measured covering factors, 
implying that it does not stemfromlow covering factors. We obtain T s = (-0.68 ±0.17) x [Z/H]+(2.13 ±0.21) 
from a linear regression analysis. Our results indicate that the high T s values found in DLAs do not arise from 
differences between the optical and radio sightlines, but are likely to reflect the underlying gas temperature 
distribution. The trend between T s and [Z/H] can be explained by the larger number of radiation pathways 
for gas cooling in galaxies with high metal abundances, resulting in a high cold gas fraction, and hence, a 
low spin temperature. Conversely, low-metallicity galaxies have fewer cooling routes, yielding a larger warm 
gas fraction and a high T s . Most DLAs at z > 1.7 have low metallicities, [Z/H] < -1, implying that the Hi 
in high-z DLAs is predominantly warm. The anti-correlation between T s and [Z/H] is consistent with the 
presence of a mass-metallicity relation in DLAs, suggested by the tight correlation between DLA metallicity 
and the kinematic widths of metal lines. Most high-z DLAs are likely to arise in galaxies with low masses 
(M vir < 10 10 5 M Q ), low metallicities ([Z/H]< -1, and low cold gas fractions. 
Subject headings: galaxies: evolution — galaxies: ISM — radio lines: galaxies 



1. INTRODUCTION 

Damped Lyman-a systems, with Hi column densities 
Nhi > 2 x 10 20 cm" 2 , are the high-z gas-rich counterparts of 
today's normal galaxies, and crucial to understanding galaxy 
evolution. However, despite much observational effort, our 
knowledge of the typical size, structure and internal condi - 
tions of high z DLAs is yet limited (e.g. I Wolfe et all 120051) . 
A controversial issue is the Hi temperature distribution in the 
absorbers, whether most of the Hi is in a warm, low-density 
phase (the warm neutral m edium, WNM, with k inetic temper- 
ature T k ~ 5000-8000 K; IWolfire et alj|1995l) . or whether a 
significant fraction is in a high density, cold phase (the CNM, 
with 7j ~ 40-200 K). For DLAs towards radio-loud quasars, 
this can be determined by combining the optical depth in the 
redshifted Hi 21cm line with the Hi column density mea- 
sured from the Lyman-a profile to obtain the column-density- 
weighted harmonic mean spin temp erature T s along the sight- 
line (e.g. iKanekar & Briggs 2004). High-z DLAs have been 
show n to have higher T s values (> 700 K; e.g. IWolfe & Davis! 
Il979t ICarilli et all [1996; Kanekar & Chengalur 2003]) than 
those ty pically found in the Milky Way or local spirals (T s < 
300 K: iBraun & Walterbosl 119921) . The simplest interpreta- 
tion of this, in the context of stable two-phase models, is that 
high-z DLAs conta in larger WNM fractio ns than local disks 
(Carill i et al|1996l:IChengalur & Kanekarll2000l) . Conversely, 
I Wolfe et all (120031) argue that the detection of CII* absorption 
in a number of high-z DLAs indicates sizeable CNM fractions 
(~ 50%) in half the DLA population. 

A plausible cause for the putative higher WNM fractions in 
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high-z DLAs is their typically-low metallicity [Z/H], imply- 
ing fewer routes for gas cooling. If high T s values in DLAs 
arise due to low absorber metallicities, one would expect an 
anti-correlation between 7^ and [Z/H], with low- 7/^ DLAs hav- 
ing h igh metallicities, and vice-versa (IKanekar & Che ngalur 
2001). We report here the detection of the predicted anti- 
correlation between T s and [Z/H], supporting the conclusion 
that Hi in high-z DLAs is predominantly in the warm neutral 
medium. 

2. THE SAMPLE 

Over the last decade, we have carried out Hi 21cm absorp- 
tion studies of DLAs towards compact , radio-loud quasars to 
measure their spin temperat u res (e . g. IKanekar & Che ngalur 
12001 IKanekar et al.1 120061 l2l)07t lYorket all 120071) . and 
have also measured the DLA covering factors through low- 
frequency very long baseline interferometry (VLBI) studies 
( Kane kar et alj|2009al) . The VLBI images yield the fraction 
of compact radio emission, and thus a lower limit to the DLA 
covering factor. However, no additional radio emission is de- 
tected up to scales of ~ 1", indicating that the remaining emis- 
sion arises from much larger scales (> 10 kpc), and is unlikely 
to be covered by the foreground DLA. The radio core fraction 
thus provides a good estimate of the DLA covering factor. 

We have also obtained metallicity estimates for most of the 
DLAs with Hi 21cm studies from our own observations or the 
literature. There are 26 DLAs, at 0.09 < z < 3.45, with esti- 
mates of both T s and [Z/H], of which 20 have covering factor 
estimates from low-frequency VLBI studies, and 21 have esti- 
mates of dust depletion, [Z/Fe]. The 26 absorbers of the sam- 
ple are listed in TableQ] whose columns contain (1) the quasar 
name, (2) the DLA redshift, (3) the Hi column density mea- 
sured from the Lyman-a profile, (4) the DLA covering factor 
/, (5) the spin temperature T s , or, for Hi 21cm non-detections, 
the 3<T lower limit to T s (taking into account the DLA cov- 
ering factor, when known), (6) the metallicity [Z/H], (7) the 
dust depletion factor [Z/Fe], (8) the transitions used for the 
[Z/H] and [Z/Fe] estimates, and (9) references for [Z/H] and 
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TABLE 1 
The full sample of 26 DLAs. 



QSO 




Nhi 
x 10 20 cm -2 


f 


T s 
K 


[Z/H] 


[Z/Fe] 


Z,Fe 


Ref. 


0738+313 


0.0912 


15±2 


0.98 


775 ±100 


< -1.14 


< 0.48 


ZnFe 


1 


0738+313 


0.2212 


7.9± 1.4 


0.98 


870 ± 160 


< -0.7 


< 0.74 


Zn,Cr 


1 


0952+179 


0.2378 


21.0±2.5 


0.66 


6470 + 965 


< -1.02 


< 0.63 


Zn,Cr 


1 


1127-145 


0.3127 


51 ±9 


0.9 


820 ± 145 


-0.90 ±0.11 




Zn,- 


2 


1229-021 


0.3950 


5.6 ± 0.9 


0.42 


95 ± 15 


-0 45 ± 14 


> 0.83 




3 


0235+164 


0.5242 


50 ± 10 


1.0 


210±45 


-0.14±0.17 


1.73 ±0.24 


X-ray,Fe 


4 


0827+243 


0.5247 


2.0 ±0.2 


0.7 


330 ± 65 


-0.62 ±0.05 




Fe+0.4,- 


1 


1122-168 


0.6819 


2.8 ± 1 




> 1445 


< -1.47 


< -0.15 




1,2 


1331+305 


0.6922 


17.8±0.8 


0.9 


965 ±105 


-1.35±0.03 


0.28 ±0.03 


Zn Fe 


5 


0454+039 


0.8596 


4.9 ±0.2 


0.5 


> 690 


— 99 ± 12 


0.00 ±0.14 




1 


2149+212 


0.9115 


5 ± 1 




> 2700 


< -0.93 




Zn- 


6 


1331+170 


1.7764 


15 ± 1.4 


0.72 


1260 ±335 


-1.20 ±0.04 


0.847 ±0.008 




7 


1157+014 


1.9436 


63 ± 15 


0.63 


1015 ±255 


-1 40 ±0 10 


38±0 16 


ZnFe 


g 


0458-020 


2.0395 


60 ± 10 


1.0 


560 ±95 


-1.27 ±0.07 


0.47 ± 0.05 


ZnFe 


7 


0311+430 


2.2898 


2.0 ±0.5 




140 ±35 


> -0.6 




Si.Fe 


9 


0432-440 


2.3021 


6.0 ± 1.4 




> 995 


-1.12±0.15 


33±0 21 


Si,Fe 


10 


0438-436 


2.3474 


6.0 ± 1.4 


0.59 


900 ±250 


-0.68 ±0.15 


0.62 ±0.21 


Zn,Fe 


10 


0405-331 


2.5693 


4.0 ±0.9 


0.44 


> 785 


-1.40±0.15 


0.34 ±0.21 


Si,Fe 


10 


0913+003 


2.7434 


5.5 ± 1.3 




> 800 


-1.47 ±0.15 


0.12±0.21 


Si,Fe 


10 


1354-170 


2.7799 


2.0 ±0.7 




> 795 


-1.86±0.16 


0.67 ±0.10 


Si,Fe 


7 


2342+342 


2.9084 


13 ±3 


0.71 


> 1705 


-1.23±0.15 


0.36 ±0.1 3 


Zn,Fe 


11 


0537-286 


2.9742 


2.0 ±0.5 


0.47 


> 520 


< -0.44 




Zn,Fe 


10 


0336-017 


3.0619 


15.8±3.6 


0.68 


> 6670 


-1.40 ±0.10 


0.37 ±0.03 


S,Fe 


7 


0335-122 


3.1799 


6.0 ± 1.4 


0.62 


> 1850 


-2.56±0.15 


0.05 ±0.21 


Si,Fe 


10 


0201+113 


3.3875 


18±3 


0.76 


1050 ±175 


-1.26±0.15 


0.18±0.21 


S,Fe 


12 


1418-064 


3.4482 


2.5 ±0.6 


0.69 


>680 


-1.48±0.15 


0.24 ±0.21 


Si,Fe 


10 



c ontain the or iginal refer ences; 

(3) IBoisse et all 119931 : (4) Junkkarinen et al. (2004); (5) Wolfe et al. (2008); (6) Nestor et al. (2008); (7) Prochaska et al. (2007); (8) ILedoux et~af] J2006T) ; 



all values are scaled to the solar abundances of Lodders (2003). (1) Kulkarni et al. 12005); (2) Kane kar et al. (in prep.); 

etaUCTOl); ( 8)^" 

(9) Ellison et al. (2008); (10) Akerman et al. (2005); (11) Prochaska et al. (2003b); (12) Ellison et al. (2001). 



[Zn/Fe] values. In all but two cases, Z = Zn, S, Si, in order of 
preference. The exceptions are the z ~ 0.524 DLAs towards 
0235+164 and 0827+243, wher e the metallicities are, res pec- 
tively, from an X-ray spectrum (Junkkarinen et al. 2004]) and 
[Z/H]=[Fe/H]+0.4 [following IProchaska et al.l (l2003al) l. The 
sample contains 14 T s measurements and 12 lower limits, and 
20 [Z/H] measureme nts, five upper limit s, and one lower limit 
(towards 0311+430; Ellis on et al]|2008l) . Detailed references 
for the metallicities and spin temperatures are provided in 
Kanekar et al. (in prep.). 

Note that the sample does not include four systems where 
the radio emission is clearly extended on scales >> 1", and 
the covering factor is likely to be low, / < < 1 ; these are the 
Hi 21 cm absorbers at z ~ 0.437 towards 3C196 (Brigg set al.l 
12001b and z ~ 0.656 towards 3C336 (ICurran et al l 12007a) 
(where the Hi 21cm absorption arises towards extended 
lobes, with little radio flux density associated with the quasar 
core), and the Hi 21cm non-detections at z ~ 1.3911 to- 
wards QSO 0957+561 and z ~ 1 .4205 towards PKS 1354+258 
dKanekar & Che ngalur 2003). We have also excluded "asso- 
ciated" DLAs, lying within ~ 3000 km s -1 of the quasar, as 
conditions in these absorbers could be affected by their prox- 
imity to an active galactic nucleus. 

Preliminary results of this study showing the anti- 
correlation between T s and \ZfH] were pre sented in 
iKanekar & Chengalur! (120051) . ICurran et all (l2007bl) later also 
found weak evidence for the anti-correlation, but did not have 
covering factor estimates for most high-z DLAs, and hence 



could not rule out the effects of unknown covering factors. 
The results of the present Letter are based on new Hi 21cm 
data on 1 1 absorbers, and VLBI estimates of the DLA cover- 
ing factor for most systems, and yield the first clear evidence 
for a relation between the metallicity and the Hi temperature 
distribution in the absorbers. 

3. RESULTS: AN ANTI-CORRELATION BETWEEN 7} AND [Z/H] 

Fig. GlA] shows T s plotted versus [Z/H] for the 26 DLAs 
of the full sample; it is apparent that low T s values (< 300 K) 
are obtained at high metallicities ([Z/H]> -0.6), while high 
T s values > 700 K are found at low metallicities, [Z/H]< -1. 
We used the non-parametric generalized Kendall ra nk correla- 
tion statistic dBrown et al.ll 19741 llsobe et al.l fl986). as imple- 
mented in the ASURV package (the BHK statistic), to test for 
a correlation between T s and [Z/H], treating the latter as the in- 
dependent variable. This allows limits in both variables to be 
treated consistently. Errors on individual measurements were 
handled through a Monte-Carlo approach, using the measured 
values and la errors for each absorber to generate 10 4 sets of 
26 pairs of [Z/H] and T s values. The statistical significance 
of the BHK statistic quoted below is the average of values 
obtained in these 10 4 trials. For the full sample, the anti- 
correlation is detected at ~ 3.6<r significance; the probabil- 
ity of this arising by chance is ~ 3 x 10~ 4 . If we exclude the 
z = 0.524 DLA towards 0235+164 (whose metallicity is from 
an X-ray study), and use [Fe/H] as a lower limit to [Z/H] for 
the DLA towards 0827+243, the anti-correlation is detected 
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FIG. 1. — The spin temperature T s plotted against metallicity, [Z/H], for [A] (left panel) the 26 DLAs of the full sample, and [B] (right panel) the ten DLAs 
with measurements of both T s and [Z/H]. The dashed line shows the linear fit to the relation between [Z/H] and T s . 



at ~ 3.1a significance, showing that it does not arise from 
incorrect metallicity estimates in these absorbers. 

A linear regression analysis was used to obtain the best- 
fit relation between T s and [Z/H], using the ten systems 
with measurements of both T s and [Z/H] (i.e. excluding 
all limits), all of which also have covering factor estimates; 
these are shown in Fig. TfBI. The BCES(Y/X) estimator 
dAkrita s & Bershady Il996l) was used for this purpose, treat- 
ing [Z/H] as the independent variable X; this method takes 
into account measurement errors on both variables, as well 
as the possibility that these errors are correlated (which ap- 
plies here, as 7^ and [Z/H] are both derived from Nht). We 
obtain T s = (-0.68 ± 0.17) x [Z/H] + (2.13 ± 0.21); the fit 
is shown as a dashed line in Figs. [T|A] and [B]. Exclud- 
ing the DLAs towards 0235+164 and 0827+243 yields 7; = 
(-0.77 ±0.40) x [Z/H]+(2.02 ±0.48), from eight DLAs with 
measurements of T s and [Z/H]. 

We also tested for a correlation between the dust depletion 
factor [Z/Fe] and T s , using the 21 systems with estimates of 
both quantities. The anti-correlation between [Z/Fe] and T s 
has ~ 2.2cr significance, reducing to ~ 1.6cr significance on 
excluding the possibly-unreliable X-ray metallicity estimate 
in the DLA towards 0235+164. We thus do not find significant 
evidence for a relation between [Z/Fe] and T s , although this 
s ub-sample contains fe w DLAs with low spin temperatures. 

ICurran et al.l (12007b!) note that an observed anti-correlation 
between T s and metallicity might arise due to unknown low 
covering factors. We hence estimated the BHK statistic for 
the sub-sample of 20 DLAs with covering factor estimates 
dKanekar et al.ll2009al) . The anti-correlation between T s and 
[Z/H] is then detected at ~ 3a significance, with a probability 
of ~ 0.003 of chance occurrence. Low covering factors are 
thus unlikely to be the cause of the anti-correlation between 
T s and [Z/H]. 

It is possible that the observed anti-correlation between T s 
and [Z/H] is not a "primary" relation, but stems from an un- 
derlying relation between other physical quantities (e.g. ve- 



locity spread and metallicit y; iLedoux et al.1 120061) . For ex- 
ample, ICurran et al.l (12007b!) argue that an observed correla- 
tion between T s and [Z/H] might arise due to a relation be- 
tween the velocity spread of Hi 21cm absorption AV21 and 
the rest equivalent width of the MgllA2796 line, Wq 2796 , 
each tracing the dynamics of the absorber. They note that 
T s oc N H i/ [/ r 2I dV] w N ffl / [ 

7"2i,maxAV2i] > suggesting that 
T s oc I/AV21. W A2796 is known to correlate with metallicity 

dMurphy et al.ll2007l) . so, if AV 2 i also correlates with W^ 2796 , 
the three relations could yield an observed anti-correlation be- 
tween T s and [Z/H]. However, the Hi column density in DLAs 
does not correlate with either [Z/H] or W A2796 (e.g. iRao et all 
2006), and the DLAs of Table Q] have Nhi values extending 
over « 1 .5 dex. If the T S -[Z/H] relation arises due to an un- 
derlying relation between AV21 and W A2796 , this spread in 
Nhi values would imply that the T S -[Z/H] relation would be 
weaker than the relation between AV21 an d VV A2796 (and vice- 
versa, if the T S -[Z/H] relation is the "primary" relation). In 
other words, the "primary" relation would have the highest 
statistical significance in samples of similar size, as uncor- 
rected (and variable) parameters like Nhi would dilute the 
si gnificance of any deriv ed relations. 

Kane kar et al.l d2009bl) used the BHK test to test for corre- 
lations between W A27% and both the Hi 21cm velocity spread 
AV21 and the integrated Hi 21cm optical depth J T2idV, in 
a sample of Mgll absorbers and DLAs. They found the pu- 
tative correlation between AV21 and W A2796 to have ~ 1 J a 
significance (with 23 absorbers), and that between J T2\dV 
and Wq 2796 to have ~ 2.2a significance (with 24 absorbers). 
Both these trends are significantly weaker than the relation be- 
tween T s and [Z/H] found here (~ 3.6a significance), in sam- 
ples of similar size. It is thus unlikely that the anti-correlation 
between T s and [Z/H] arises from an underlying relation be- 
tween Wq 2796 and either AV21 or J T2idV. The present data 
indicate that the anti-correlation between T s and [Z/H] is the 
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"primary" relation, given its higher statistical significance. 

4. DISCUSSION 

High spin temperature estimates in high-z DLAs have 
been the focus of much debate, with suggestions that these 
might arise due to different sightlines at optical and ra- 
dio wavelengths ([Wolfe et al.1 120031) or low covering factors 
dCurran et alj|2005l) . It is possible that individual T s measure- 
ments might be in error due to these effects. However, recent 
low-frequency VLBI studies have shown that high T s values 
are not the result of low DLA covering factors (Kanek ar et alJ 
l2009al) . Further, the anti-correlation between T s and [Z/H] 
is detected in the sub-sample of DLAs with covering factor 
measurements, indicating that this relation too is not caused 
by low covering factors. Finally, differing optical and radio 
sightlines (e.g. if the radio emission is extended on scales 
larger than the size of the optical quasar) would result in T s 
estimates that are uncorrected with the "true" T s values along 
the optical sightline. This would weaken any underlying cor- 
relation between T s and a quantity like [Z/H], measured along 
the pencil beam towards the optical quasar. The detection of 
the predicted anti-correlation between T s and [Z/H] thus in- 
dicates that line-of-sight issues are also not the source of the 
observed high T s values. We conclude that the high T s esti- 
mates in high-z DLAs are most likely to arise due to larger 
WNM fractions in DLAs than typical in local spiral galaxies. 

In the local universe, the metallici ty and mass of galax - 
ies are known to be correlated (e.g. iTremonti et al] I2004I) . 
A similar correlation, between metallicity and stellar mass, 
has been found in emission-selected, high-z galaxies (e.g. 
Savag iio"et al.l 12005). E vidence for a mass- metallicity rela- 
tion in DLAs is unclear. Led oux et al.l d2006l) argue that the 
observed correlation between the kinematic width of unsat- 
urated low-ionization metal profiles AV90 and metallicity in 
high-z DLAs arises from an und erlying mass-metallicity re- 
lation. Conversely, IZwaan et al.1 (|2008) note that AV90 is a 
weak indicator of galaxy mass at z = 0, as the kinematic width 
measured along a pencil beam depends on the galaxy incli- 
nation and the impact parameter. However, the above cor- 
relation betwee n AV90 and m etallicity has been reproduced 
in simulations dPontzen et al.l 120081) . which s uggest that it 
does s tem from a mass-metallicity relation. Pontz en et al.l 
(2008) find that sightlines through DLAs with low virial 
masses, M wu < 10 9 ' 5 M Q , typically yield low metallicities, 
[Z/H]< -1.7, while the typical metallicities are [Z/H]^ -1.2 
in intermediate-mass galaxies (M vil - - io 9 - 5 " 10 - 5 M Q ), and 
[Z/H]> -1 in high-mass systems (M vir > 10 ia5 M Q ; for com- 
parison, M vir ~ 10 12 M(7) for the Milky Way ; e.g. IXue et all 
120081) . Interestingly, iProchaska et al.1 (12008) find a tight cor- 
relation between the rest equivalent width in the saturated 
SillA1526 line \V\526 and metallicity, and argue that ^1526 
tracks dynamical motions in the halos of DLAs, with the 
correlation between W1526 and [Z/H] arising due to a mass- 
metallicity relation of the same slope as that seen in low-z 
galaxies. 

If a mass-metallicity relation is present in DLAs, low- 
mass DLAs would have low metallicities. Low-mass galax- 
ies also have low thermal pressures, while a high pres- 
sure is needed for the formation of a stable mult i -phase 
medium at a low metallicity. Fig. 6 of IWolfire et al.l (11995b 
shows that, at low pressures and metallicities, Hi exists 
mostly in the warm phase. Sightlines through such galax- 
ies would thus typically yield a high spin temperature. Con- 



versely, the mass-metallicity relation implies that high-mass 
DLAs would have high metallicities, along with high ther- 
mal pressures. Such galaxies would have significant CNM 
fractions, with most sightlines yielding low spin tempera- 
tures. These effects would yield the anti-correlation be- 
tween T s and [Z/H] detected here, due to the paucity of 
cooling routes in low-metallicity galaxies (|Norman &"S paans 
[l997t IKanekar & Chengalurfl200lT; IKanekar & Briggsll2004l) . 
The anti-correlation between T s and [Z/H] is thus consis- 
tent with the presence of a mass-metallicity relation in DLAs 
dLedoux et al.ll2006t IProchaska et all 120081) ; low-metallicity 
DLAs are likely to have high T s values due to their low CNM 
fractions. 

While the relation between T s and [Z/H] is consistent with 
the presence of a mass-metallicity relation in DLAs, it is also 
possible that the r s -[Z/H] relation is a local one, arising due 
to line-of-sight issues. For disk galaxies, the cross-section 
for DLA incidence is largest for sightlines through the out- 
skirts of the galaxy. If high-z DLAs have steep metallic- 
ity gradients, such sight lines would ty pically encounter low 
metallicities (e.g. IZwaan et al.l 12005). The lack of local 
cooling routes could then result in low CNM fractions, and 
high T s values along these sightlines. Conversely, sightlines 
through the central regions of large disk galaxies would typ- 
ically have high metallicities and high CNM fractions. Such 
line-of-sight effects could thus also yield the anti-correlation 
of Fig. Q] While such steep metallicity gradients have not been 
seen in local galaxies and low-z DLAs (typical gradients are 
~ -0 03 — 0.04 dex kpc -1 ; iChen et al.l 2005; Br esolin et al.1 
2009), they are not ruled out in high-z DLAs. We note, in 
passing, that the z ~ 0.524 DLA towards 0827 +243 has one 
of the largest DLA impact parameters (~ 27 kpc; Steid el et al.l 
120021) . and yet has both a low T s (330 K) and a high metallic- 
ity, [Z/H]= -0.6. 

Fig-EA] shows that low T s values are only found in DLAs 
with [Z/H]> -0.6. Only 13 of the 1 53 DLAs at z > 1.6 in 
the sample of [Prochaska et al.l (120071) have such high metal- 
licities; the majority of high-z DLAs have [Z/H]< -1.0. The 
anti-correlation between T s and [Z/H] then implies that most 
high z DLAs have high T s values. Converse ly, six of seven- 
teen DLAs at z < 1 have [Z /H]> -0.6 (e.g. IProchaska et al.1 
l2003atlKulkarni et al.ll2005l) . While the number of metallicity 
measurements at z < 1 is yet small, the present data indicate 
a higher fraction of high-metallicity DLAs at low redshifts. If 
so, the T S -[Z/H] anti-correlation implies that the fraction of 
DLAs with high CNM fraction should also be larger at z < 1 . 
This is consistent with the observed increase i n the detection 
rate of Hi 21cm absorption in DLAs at z < 1 (Kanekar et al. 

l20Tj9bh . 

Finally. IWolfe et all (120031) found that the strong CII* lines 
detected in half their DLAs could not be explained by ab- 
sorption in pure WNM, and implied the presence of cold 
Hi. They argue that Hi in DLAs without CII* absorption is 
likely to be predominantly warm, but that a two-phase model 
with a CNM fraction of ~ 50% is consistent with the ob- 
served bolometric luminosity in systems with CII* detections. 
However, lower CNM fractions are not ruled out. For exam- 
ple, the z ~ 3.39 DLA towards PK S 0201+113 shows strong 
CII* absorption (Wolfe et al. 2003), but has a low metallicity, 



[S/H]= -1.26 (Ellison et alJpOOll) . and a low CNM fraction, 
< 17% (IKanekar et al.1120071) . This suggests that, while DLAs 
showing CII* absorption contain some CNM, the Hi content 
of low-metallicity absorbers is still likely to be dominated by 
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the WNM (see also lSrianand et al.ll2005l) . 

5. SUMMARY 

We report the detection of an anti-correlation between spin 
temperature and metallicity (with ~ 3.6a significance in the 
non-parametric BHK test) in a sample of 26 DLAs at 0.09 < 
z < 3.45. For 20 systems, the absorber covering factor has 
been estimated from low-frequency VLBI studies; the anti- 
correlation between T s and [Z/H] is detected here at ~ 3a sig- 
nificance. A linear regression analysis using the BCES esti- 
mator finds the relation T s = (-0.68 ±0.17) x [Z/H] + (2.13 ± 
0.21) between T s and [Z/H]. Low spin temperatures, T s < 
300 K, are only found in high-metallicity DLAs (with [Z/H]> 
-0.6), while high T s values (> 700 K) are obtained in low- 
metallicity ([Z/H]< -1) absorbers. The fact that a relation is 
seen between T s and [Z/H] implies that the high T s values in 
DLAs are unlikely to be an artifact arising from differences 
between radio and optical sightlines through the absorbers, 



and are instead likely to reflect the underlying gas temper- 
ature distribution. The majority of DLAs at z > 1.6 have 
[Z/H]< -1, implying that most of the Hi in DLAs must be 
in the warm phase, with small CNM fractions. The observed 
anti-correlation between T s and [Z/H] is consistent with inde- 
pendent evidence for the presence of a mass-metallicity rela- 
tion in DLAs. The majority of high-z DLAs are likely to be 
galaxies of low mass and metallicity, with most of the neutral 
gas in the warm phase. 
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